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Abstract 


Formulae  are  derived  for  the  Hough  components  and  the 
mean  rate  of  solar  radiational  heating  by  ozone  absorption 
in  the  presence  of  a  lower  reflecting  layer.  A  model  of 
ozone  densities  based  on  observational  data  is  presented 
for  each  of  the  four  seasonal  months  January,  April,  July 
and  October  at  heights  0(5)80  km  and  latitudes  0(15)90°. 
Diurnal  and  semi-diurnal  components  of  heating  are 
evaluated  for  migrating  modes  and  discussed  in  relation 
to  features  of  the  adopted  ozone  model.  Comparisons 
are  made  with  earlier  evaluations. 
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f 1.  Introduction 

i 

A  major  source  of  thermal  excitation  of  atmospheric 
tides  is  that  of  absorption  of  solar  radiation  by  ozone* 
Calculations  by  Butler  and  Small  (1963)  showed  that  heating 
in  the  region  of  50  km  altitude  contributes  significantly 
to  surface  pressure  oscillations  through  particular  modes 
of  oscillation,  notably  the  leading  semi-diurnal  (2,2,2) 
mode.  Prom  ozone  absorption  cross-sections,  solar  spectral 
data  and  height  profiles  of  ozone  concentration,  heating 
rates  may  be  calculated  for  a  given  time  of  the  year  at 
given  heights,  latitudes  and  times  of  day.  By  extracting 
the  appropriate  Fourier  time  component  (e.g.  diurnal  or 
semi-diurnal)  and  taking  the  latitudinal  average  weighted 
by  the  selected  Hough  function,  the  corresponding  component 
of  the  heating  rate  is  obtained.  Lindzen  (1967)  obtained 
diurnal  Hough  components  of  equinox  heating  from  the 
radiative-photochemical  model  of  Leovy  (1964).  Semi-diurnal 
components  have  likewise  been  derived  (Chapman  and  Lindzen, 
1970;  Lindzen  and  Hong,  1974;  Hong  and  Lindzen,  1976). 

The  purpose  of  this  paper  is  to  evaluate  the  (1,1,1), 
(2,2,2)  and  other  diurnal  and  semi-diurnal  Hough  components 
of  ozone  heating  in  each  of  the  four  seasons  of  the  year 
on  the  basis  of  observed  ozone  concentrations.  Particular 
attention  is  given  to  the  construction  of  a  seasonal- 
latitudinal  ozone  model  ($4).  For  calculating  the  solar 
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energy  absorption  at  many  heights  with  global  coverage  a 
rapid,  method  of  computation  is  required.  To  meet  this 
need  the  parametric  treatment  of  Lacis  and  Hansen  (1974-) 
has  been  adopted  and  is  described  in  .§2. 

2.  Formulation  of  heating  rate  due  to  ozone  absorption 
Let  A  be  the  ozone  absorption  at  height  z  of  solar 
radiation  incident  on  the  atmosphere  at  zenith  distance  5  , 
then  the  rate  of  heating  within  a  cylindrical  volume  of 
unit  horizontal  base  area  and  height  dz  inclined  at  an 
angle  £  to  the  vertical  is  -Sc«*£ciA  ,  where  SQ  is 
solar  energy  flux  and  dA  is  the  increase  in  A  corresponding 
to  an  increase  dz  in  z.  The  heating  rate  per  unit  mass 
of  atmosphere  is  then 

j*  te.  -  v5>  —  (2-0 

*  r 

where  is  air  density  and 

v  =  ccs  x  G*  ■  0 

SQ  is  taken  as  0.1 353  W  cm“2  (Thekaekara,  1973)  with  an 
inverse  square  dependence  on  the  Sun-Earth  distance.  By 
spherical  trigonometry 

8  S  Jbm*  8  S  Ctfi  t  Ca  •  O 


V 
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where  S  is  solar  declination,  6  is  colatitude  and  It  -ft* 
is  local  time  in  radians. 

To  evaluate  the  derivative  in  (2.1)  we  use  the 
parameterization  for  ozone  absorption  of  Lacis  and  Hansen 

(1974) 

A  -  0  0*119*.  +  i  o>ix  +  Q.Ofeygjt  C2  ^ 

°  It  O.oui^t  0.000*23*  (I  t  1 59,fe  x.)0'*0*"  I  •+  (loJ.feaO* 

where 

oc  -  u.H  CaO 

u  is  the  ozone  amount  (  in  cm  at  NTP)  in  a  vertical  column 
above  height  z  and  ft  is  the  magnification  factor  due  to 
slant  path  and  refraction 

=.  £$/ (l«4v%  I  Ca  t) 

The  first  term  on  the  right-hand  side  of  (2.4)  relates  to 
absorption  at  visual  wavelengths  and  has  a  numerical 
precision  of  nearly  four  decimals  in  the  interval 
10  ^  x  ^  10  cm.  The  second  and  third  terms  relate  to 

ultraviolet  absorption  and  are  accurate  to  within  0.5  per 
cent  in  the  interval  10  ^  x  1  cm. 

In  (2.1)  dA/dz  would  be  equal  to  dAQ/dz  in  the  absence 
of  a  reflecting  layer  beneath  the  ozone  absorbing  region. 

In  general  reflections  from  the  surface  and  lower  atmosphere 


provide  an  upward  flux  of  diffuse  radiation  and 


4A  _  JA0c»)  _  ^  4/\0C*‘)  (2.7) 

dz  d* 

where  x'  is  the  ozone  path  length  traversed  by  solar 
radiation  on  reaching  height  z  from  below.  We  have 

x'  -  ut  M  +  Cut-u)M  <?-*) 

where  u^.  is  the  total  ozone  amount  in  a  vertical  path  above 
the  reflecting  layer  (the  surface  for  clear  sky  conditions  or 
the  cloud  tops  for  cloudy  skies)  and  R  =  1,9  is  an  effective 
magnification  factor  for  the  upward  diffuse  radiation.  The 
negative  sign  arises  in  (2.7)  as  Aq(x' )  decreases  with  z 
whereas  A(x)  increases  with  z.  P.  is  the  albedo  of  the 
reflecting  layer  and  depends  on  the  effective  albedo  of 
the  lower  atmosphere  Ra  and  the  surface  albedo  R^.  It 
may  be  shown  that  (Lacis  and  Hansen,  197^) 

K  -  R«.  +  O (M> 

R_  and  hence  R  are  functions  of  V  ,  and  5_  is  defined  by 

8  cL 

_  f  I 

*  Z.  j  R^tv)  v  dv  (?  io) 

c 

Let  (1-K)  be  the  fraction  of  cloud-free  sky  having 
an  albedo  (from  above)  of  R  ■  R^  and  let  R  ■  R£  for  the 
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fraction  K  of  cloudy  sky  on  the  assumption  that  all  areas 
of  cloud  have  the  same  reflective  properties.  The  albedo 
of  the  reflecting  layer  on  neglecting  cloud  edge  effects  is 

R  -  0-*O  R,  +  (*  ») 

On  adopting  the  relations  of  Lacis  and  Hansen  (1974)  for 
the  reflectivity  of  a  cloud-free  atmosphere  Ra^  and  of  a 
cloudy  atmosphere  Rap,  we  obtain  R^,  Rp  from  (2. 9).  We  have 

Rat  -  o.2»}/ (|  +  o.srfev)  C2'2) 

Rai  —  0.13  r  /£f  +  (*•  *3) 

where  T  is  the  total  visual  optical  thickness  of  the  cloud 


layer.  By  (2.10) 

Ra,  *  0. /<♦•**  (?•»*) 

m  0. 13^ /()  -+  o.iSt)  &  I5’) 

Hence  (2.11)  may  be  written  by  (2.9)*  (2.32)  to  (2.15) 

^  *  C#/(|  +0.8/6  v)-»-  cx  (2.(0 

where 

C/»  0.215  (J- -  fc,)  -'7) 

cA—  ■+  *■  £  +  0  “  Ca.li) 

o-*./)Vcu*v**)  y* ,>a) 


cr 


>  mu',  a  U  kjfc*.  . 1 
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3 .  Relations  for  heating  rate  components 

We  express  time  (G.M.T.)  t  and  longitude  ^  in  radians 
and  expand  J  as  a  Fourier  series  in  t  in  the  interval  (0,  2ic  ) 

cl 


0*0 


w.  =1,2,..  relate  to  the  diurnal,  semi-diurnal,.,  components 
of  heating.  Defining 

Ji  -  J'J?  +  <:  xj  C3.i) 


(3.3 )  may  be  written  as 

z,  -  £  j*  eimt 

m*  0 

Since  Jm  is  periodic  in  ^  in  the  interval  (0,  2”it )  it 
may  be  expanded  as 


(.3  5) 


■rm-  E  ■**■'+) 

^  •  A  S.  'J 


0<0 


where  Jq,  j|  are  complex.  It  is  convenient  to  replace 
Jq,  j|  by  J^s  defined  by 


<  »  J-c’ 

=  J-c"  +  iJ/ 

(->  >0) 

Then  (3.4-)  becomes 

jr  =  £  j*  -L*+ 


C36) 


•  *v. 
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and  (3.3)  becomes 


7  =  Itl  E  T-* 


C*7) 


is  a  function  of  colatitude  8  which  is  expanded  for 
m  0  in  terms  of  normalized  Hough  functions  flu  C|0 


J'  *  I  J“  0"  ((*) 

w  n.  ifft*  1 


(»w-4  (?)  (.J.S) 


where 


JA  =  Crt  0 


(?-l) 


The  summation  in  (3*8)  is  taken  over  all  members  of  the  set 
of  eigenfunctions  of  Laplace's  tidal  equation  for  given 
w*  (*0)  and  s.  The  full  expansion  of  J  is  then 


It  may  be  shown  using  (3*  "  that 


*1*  iff 


*  4tta  j  j 


C3.//; 


“  air*  j  (  **  *  dt  J<f>  L*  +  °) 


and  on  using  the  orthogonality  properties  of  ^  that 
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for  wi  f  0 


7^  =  — 

m.n.  3.1* 


fY'f1"  _  a-W 

J  I  J  <(“>, t)  0_J|O  «.  JtJfJp  (?••*) 

-•  Jt»  '\,  / 


Local  time  in  radians  is  expressed  by  either  side  of  the 


equation 


(^  +  t  *  IT  +  t 


(3.  .3) 


Hence  from  (2.1),  (2.7)>  (2.16),  (3.11)  and  (5.12) 


T'*  .  C  ,  ** 

4tt  L  '  ° 


=  s. 


vn 

i  f.  •" 


0* 


(x3  ^  O  ) 


TTt  *T\ 


=  £ 


?  j  (  ^  two) 

'*  6  ’  o-'O 


where 


-i  "o 


r-*\  * ,  / 

>L  V  <2.  ctl 


C?  *7) 


H"F  -  Cf>  +  c,  -  c/.yv  (j<0 

f.<*r.-rt  =  -  f  »  J*' a’  on) 

J-K  ** 

{,(*,«.<*>)  =  j1'  l^«(5t';v  (I  +  0.»I6  v)"’ (.US) 


^  dAaU‘)  .  -wt' 


V  €.  dt 


(  3.  *1) 


and  -t  ,  t  are  the  values  of  t'  at  sunrise  and  sunset 
o  o 


11 


respectively.  If  no  daily  variation  in  ozone  concentration 
is  present,  u  in  (2.5)  is  independent  of  t1  and  t*-  may 

,  i 

be  replaced  by  c«ml  as  x,  x1  and  V  are  even  functions 


of  t‘. 


If  ozone  and  air  density  are  independent  of  longitude, 


(3.14)  and  (3.35)  may  be  written  as 

j.-v  -  is.?: 

- 


(4=0) 
(**  ■=♦=  0) 


(l.iu) 


t»+0) 


where 


s:  -  c  +  (j.«) 


I 


0 

i  LC'<s-*n)A 

*«Lcj€'  Jt 


C3.23) 


y  *'.*) 


4.  A  latitudinal  and  seasonal  ozone  model 

Table  1  presents  a  model  of  ozone  concentrations  that 
has  been  prepared  from  observational  data.  At  high  latitudes 
and  at  heights  above  50  km  observational  support  is  limited 
and  values  are  of  a  provisional  nature.  Belcw  50  km,  data 
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Table  1.  (continued)  April  N  (or  October  S)  and  October  N  (or  April  S) 
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were  sufficiently  abundant  for  only  those  dates  falling 
within  the  months  of  January,  April,  July  and  October  to  be 
used,  but  at  greater  heights  data  for  dates  up  to  one  month 
on  either  side  of  these  months  were  included. 

A  longitudinal  ozone  dependence  has  been  apparent  for 
some  time  from  N  hemisphere  total  ozone  amounts  (London, 

1963)*  The  main  feature  of  .the  distribution  is  the 
existence  of  three  high-latitude  regions  of  maximum  ozone 
over  N.E.  America,  E.  Asia  and  N,  Europe  which  become 
particularly  well  developed  and  steady  in  late  winter 
and  early  spring  (Khrugian,  1973) •  Over  the  last  10  years, 
satellite  techniques  have  contributed  increasingly  to  ozone 
measurement  (Krueger  et  al.,  1980).  Global  distributions 
of  total  ozone  (Prabhakara  et  al.,  1976;  Hilsenrath  et  al., 
1979)  show  good  agreement  with  ground-based  measurements 
and  have  improved  geographical  coverage  particularly  in 
the  S  hemisphere.  Evidence  of  longitudinal  variations  in 
stratospheric  ozone  has  also  been  provided  by  satellite 
observations  (Krueger  et  al.,  1973;  Frederick  et  al.,  1977); 
and  a  maximum  variation  at  35  km  of  about  12  per  cent  has  been 
observed  by  the  0G0-4  satellite  (London  et  al.,  1977)* 

Such  data  are  however  inadequate  for  extending  the  present 
model  in  longitude. 

The  extended  network  of  ground  stations  operating 
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since  the  I.G.Y.  has  shown  considerable  seasonal  asymmetries 
in  total  ozone  amounts  between  the  two  hemispheres  (Duetsch, 
1974).  Satellite  measurements  have  generally  confirmed 
these  results  with  improvements  in  time  resolution  and 
S  hemisphere  coverage  (Ghazi  et  al.,  1976;  Hilsenrath  et  al., 
1979).  For  vertical  ozone  distributions  however  insufficient 
data  are  available  for  a  comparable  analysis  of  seasonal 
asymmetries.  The  respective  summer,  autumn  and  winter 
distributions  have  therefore  been  represented  by  the  same 
model  for  either  hemisphere:  for  spring,  the  model  takes 
account  of  the  higher  values  at  high  N  latitudes  than  at 
corresponding  S  latitudes  with  entries  designated  A  (April  N) 
and  0  (October  S)  in  Table  1. 

The  model  between  60  N  and  60  S  is  that  prepared  by 
Lucas  (1978)  apart  from  minor  modifications  at  60°  latitudes 
and  the  addition  of  heights  0,  5  and  80  km.  Up  to  45  km, 
the  model  is  based  primarily  on  Umkehr  measurements  and 
compares  well  with  earlier  analyses.  At  35  and  40  km  the 
model  depends  equally  on  satellite  data  and  Umkehr  observations 
Above  55  km,  data  are  extremely  sparse  and  the  latitudinal 
dependence  was  obtained  by  fitting  a  lowest  order  polynomial, 
which  was  often  a  straight  line.  A  more  detailed  account 
of  the  sources  of  data  and  the  preparation  of  the  model 
for  latitudes  60  N  to  60  S  has  already  been  given  (Lucas,  1978) 
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For  extending  the  model  to  90°  latitude,  the  following 
sources  have  provided  data:  Ramanathan  and  Kulkarni  (I960); 
Krueger  et  al.,  (1975);  Duetsch  (1974);  London  et  al.,  (1977); 
and  Wilcox  et  al.,  (1977)*  Between  60  and  80  km  the  entries 
in  Table  1  are  a  simple  extrapolation  of  the  values  at  lower 
latitudes  followed  by  smoothing  in  height.  Below  90  km 
just  over  half  the  entries  in  Table  1  are  without  observational 
support  during  the  season  to  which  they  apply. 

Values  from  Table  1  at  49°  latitude  are  compared  in  Table 
2  with  those  of  the  mid-latitude  animal  mean  model  that  has 
been  adopted  for  the  1976  U.S.  Standard  Atmosphere  (Krueger 
and  Minzner,  1976).  On  averaging  the  January,  April,  July 
and  October  values  of  the  present  model  at  49°  latitude  a 
profile  is  obtained  which  agrees  with  the  standard  to  within 
2  pei"  cent  at  40  to  70  km.  In  view  of  the  appreciable 
seasonal  variation  at  the  greater  heights,  such  agreement 
is  remarkably  good.  Between  19  and  29  km  the  standard 
values  exceed  the  average  of  the  seasonal  values  by 
amounts  of  up  to  JO  per  cent  (reached  at  20  km).  At 
lower  heights  the  average  exceeds  the  standard  having  a 
sharp  increase  in  ozone  concentration  from  4  to  8  km, 
whereas  a  similar  increase  occurs  in  the  standard  values 
between  8  and  12  km:  consequently  values  differ  considerably 
at  10  km  (Table  2).  Such  differences  are  however  largely 


Ht 

km 

Jan 

at 

Apr  Jul 

45  N  latitude 

Oct 

Average  of 

4  seasons 

U.S.  Standard 
Atmosphere  1976 

10 

24.0 

27.5 

18.0 

14.0 

20.9 

11.3 

x  1011 

20 

35.5 

45.0 

36.5 

32.5 

37.4 

47.7 

x  1011 

30 

26.5 

29.5 

28.0 

28.5 

28.1 

25.2 

x  1011 

40 

61.7 

63.7 

54.2 

66.1 

61.4 

60.7 

x  1010 

50 

80.2 

63.7 

52.0 

65.0 

65.2 

66.4 

x  109 

60 

103.5 

79.4 

50.4 

56.8 

72.5 

73.3 

x  108 

70 

78.3 

68.1 

25.3 

49.5 

55.3 

54 

X 

*— 

O 

Total 

9.48 

10.89 

8.60 

8.24 

9.30 

9.27  x 

l018cm"2 

ozone 

0.353 

0.405 

0.320 

0.307  0.346 

0.345 

atm-cm 

self-balancing  in  a  vertical  integration  of  the  ozone 
concentration  as  total  ozone  amounts  agree  very  closely 
(Table  2). 


Variations  of  ozone  concentration  with  time  of  day 
are  shown  by  the  photochemical  analysis  of  Park  and  London 
(1974-)  to  be  insignificant  below  50  km,  while  above  50  km 
they  increase  with  height  and  would  appear  likely  to  effect 
the  present  calculations  at  60-65  km  and  above.  A  procedure 
has  therefore  been  adopted  whereby  the  ozone  concentrations 
of  Table  1  at  50  to  80  km  are  multiplied  by  a  height -dependent 
factor  representing  the  relative  ozone  variation  between 
sunrise  and  sunset.  The  factor  has  been  based  on  the 
calculations  of  Park  and  London  (1974- )  for  July  30  N  and 
is  applied  to  other  seasons  and  latitudes  with  rough 
approximation  by  expressing  its  time-dependence  as  a 
fraction  of  the  sunrise-sunset  interval  instead  of  local  time. 

The  observationally  derived  values  of  Table  1,  being  based 
on  techniques  limited  to  daylight  conditions,  are  taken  to 
correspond  to  the  daytime  maximum  of  ozone  concentration 
shown  by  the  theoretical  results,  i.e.  the  factor  is 
normalized  at  each  height  to  a  maximum  value  of  unity 
(Lucas,  1978). 

A  seasonal-latitudinal  model  f or  air  density  , 

which  appears  in  (3.36)  and  (3.22),  has  been  employed 
by  taking  values  from  CIRA  (1972)  between  25  and  80  km 
altitude  and  from  the  U.S.  Standard  Atmosphere 
Supplements,  1966  (COESA,  1967)  at  lower  heights. 

\ 
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5.  Surface  and  lower  atmosphere  reflectivities 

A  simple  representation  has  been  devised  for  surface 
albedo  R  by  assigning  one  of  the  values  0.07,  0.14,  0.21 

C 

or  0.76  to  each  10°  latitude  x  10°  longitude  area  according 
to  whether  the  area  is  categori'/.ed  as  all  sea,  half  sea  and 
half  land,  all  land  or  snow  covered.  No  account  is  taken 
of  daily  or  seasonal  variations  of  albedo.  The  albedo  of 
a  water  surface  depends  strongly  on  solar  zenith  angle  and 
the  value  0.07  is  that  adopted  by  Katayama  (1966)  for  diffuse 
solar  radiation:  the  same  value  was  taken  by  Lacis  and  Hansen 
(19749 •  Values  of  0.05  and  0.16  for  the  average  albedos 
of  sea  and  land  surfaces  respectively  have  been  assumed  by 
Yamamoto  et  al.,  (1974);  and  surface  albedos  of  0.1  and  0.4 
have  been  taken  to  be  representative  of  ocean  surface  and 
bright  sand  (or  old  dirty  snow)  by  Lion  et  al.,  (1978). 

Land  albedos  depend  on  the  nature  of  the  surface  and  vary 
seasonally  when  vegetated  or  subject  to  snow.  A  land 
albedo  of  0.21  has  been  adopted  here  as  this  value  is  an 
average  of  the  values  given  by  Katayama  (1966)  for  different 
types  of  land  surface.  The  value  of  0.76  for  snow-covered 
areas  is  the  mean  surface  albedo  calculated  by  Katayama 
(1967)  for  60  N  in  January. 

For  simplicity  (2.11)  was  formulated  cn  the  assumption 
that  all  cloudy  areas  have  the  same  reflective  properties. 


Katayama  (1967)  has  calculated  an  average  albedo  for  cloudy 
areas  of  0.60  (0.62  for  January  and  0.58  for  July)  which 


corresponds  to  the  quantity  R above  being  the  sum  of  cloud 
albedo  and  a  small  contribution  due  to  atmospheric  scattering 
above  the  cloud  layer  (amounting  to  about  0.045).  By 
(2.15) »  the  corresponding  value  of  optical  thickness  is  11.5. 

For  cloud  cover,  the  four  seasonal  maps  of  average 
daytime  cloudiness  prepared  by  Clapp  (1964)  from  Tiros  IV, 

V  and  VI  nephanalyses  were  adopted.  The  maps  cover  60  N 
to  60  S  and  all  longitudes,  cloud  amounts  being  shown  as 
either  <  0.50,  0.50-0.7!?  or  >  0.75.  To  each  10°  latitude 
x  10°  longitude  area  one  of  the  three  cloudiness  values  0.25, 
0.625  or  0.875  was  allocated  on  the  basis  of  these  maps.  At 
70-90°  latitudes  values  were  assigned  on  the  basis  of  the 
maps  of  Shaw  (1942). 

6.  Ozone  heating  rates 
6.1  Mean  heating 

As  a  preliminary  calculation,  mean  heating  rates  at  0° 
latitude  have  been  determined  for  the  April  ozone  profile  and 
different  solar  zenith  angles.  Results  are  plotted  in  Fig.  1 
of  Ja/cp  in  units  of  °/day  where  cp  is  the  specific  heat  of 
air  at  constant  pressure  and  J  is  obtained  from  (2.]). 

Results  are  shown  for  two  types  of  lower  reflecting  layer: 

(i)  clear  sky,  reflexion  arising  jointly  from  atmosp'heric 


Height  profiles  of  equatorial  heating  rates 
for  April  ozone  values.  Solar  zenith 
angles  are  marked,  on  curves.  Key:  -  -  -  - 
clear  sky;  -  cloudy  sky. 
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scattering  and  the  Earth's  surface  for  which  an  albedo 
of  0.14  is  taken. 

(ii)  100  per  cent  cloud  with  an  albedo  of  0.60  which  in 
conjunction  with  the  Earth's  surface  albedo  of  0.14 
gives  a  total  albedo  of  0.62  by  (2.9). 

For  an  overhead  Sun,  a  maximum  heating  rate  of  nearly 
23°C/day  is  obtained  at  44  km  (Fig.  i),  which  is  in  close 
agreement  with  previous  calculations  for  equatorial 
conditions  (Lacis  and  Hansen,  1974).  The  presence  of 
clouds  has  a  relatively  large  effect  at  lower  heights  (Fig.  l), 
e.g.  at  20  km  reflexiox:  from  clouds  increases  the  heating 
rate  by  46  per  cent  and  at  30  km  the  increase  is  20  per  cent; 
but  at  greater  heights  heating  by  incoming  UV  radiation 
dominates  that  by  reflected  radiation  which  has  traversed 
longer  absorption  paths  and  is  at  visual  wavelengths. 

Mean  heating  rates  J°/Cp  have  been  calculated  by  (3.20) 
for  0°  solar  declination  (and  April  ozone  and  air  densities) 
and  are  shown  in  Fig.  2  at  various  N  latitudes  for  case  (i) 
of  clear  skies.  From  Fig.  2  it  is  noted  that  the  highest 
heating  rate  occurs  at  30°  N  latitude  and  that  at  heights 
close  to  30  km  even  60°  N  heating  rates  may  exceed  equatorial 
values.  Such  results  would  eppear  to  be  a  consequence  of 
the  decrease  in  ozone  concentrations  towards  low  latitudes 
that  is  found  at  40-30  km  in  April  (Table  1).  The 


1  ozone  and 


observational  evidence  at  45  km  for  this  decrease  is  shown 
in  Fig.  3,  where  the  continuous  line  corresponds  to  the  values 
of  Table  1.  The  low  latitude  ozone  minimum  is  found  in 
back-scatter  UV  spectrometer  observations  from  the  Nimbus  4 
and  0G0-4  satellites  and  has  some  support  from  the  rather 
scattered  assemblage  of  Umkehr  results. 

6.2  Diurnal  Houghx -components 

\ 

Diurnal  Hough  components  of  the  heating  rate  have  been 
calculated  from  (3.21)  for  the  ozone  and  air  density  models 
described  in  §  4  and  the  reflectivity  data  of  §5. 

Normalized  Hough  functions  are  obtained  as  solutions  of 
Laplace's  tidal  equation,  methods  of  solving  which  have 
previously  been  reviewed  (Groves,  1979).  Diurnal  Hough 
functions  for  s  =  1  have  previously  been  presented  graphically 
(Groves,  1975),  but  may  differ  in  sign  from  those  used  here 
which  are  positive  at  the  equator  if  equatorially  symmetric 
or  increase  at  the  equator  with  latitude  if  equatorially 
asymmetric . 

With  s  =  m  =  1,  it  follows  from  (3*25)  that  CQ,  C^, 

are  real;  and  hence  by  (3*22)  that  is  real  provided 

x  ,n 

fQ,  f2  are  real*  (5.17)  t0  (5.19),  fQ,  f-p  f2  are 

real  if  the  daily  variation  in  ozone  concentration  is 
symmetric  v/ith  respect  to  local  noon,  but  as  described  in 
$  4  the  daily  ozone  variation  is  in  general  asymmetric  with 
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April  October 

Fig.  3  Ozone  concentrations  at  45  km  altitude  for 
April  and  October  from  0  to  60  N  latitude. 
Key:  o  Nimbus  4  BUV  spectrometer  (Krueger 
et  al.,  1973);  x  Areas  optical  ozonosonde 
(Krueger, 1973 ) ;  •  0G0-4  BUV  spectrometer 

(Anderson  et  al.,  1969); - OGO-4  BUV 

spectrometer  (London  et  al.,  1977);  ■■  -  - 
model  values  (Table  1).  Vertical  lines 
refer  to  Umkehr  determinations  and  extend 
to  one  standard  deviation  from  the  mean. 
Seasonal  symmetry  is  assumed  between  N  and 
S  hemispheres  (Lucas,  1978). 
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reszecx,  to  local  noon  above  50  km  and  J?"  then  has  an 

i ,  n 

imaginary  part.  When  calculated,  however,  imaginary  parts 
are  very  small  (of  order  10”^)  relative  to  real  parts  and  will 
therefore  be  neglected. 

The  results  obtained  for  the  three  symmetric  Hough 
components  of  heating  J,  .  ,  J?"  ,,  J?"  ^  are  shown  in  Fig.  4, 
and  those  obtained  for  _2»  Jp  ,  J^  which  are  symmetric 

components  that  are  strongly  affected  by  high-latitude 
heating  are  shown  in  Fig.  5»  Values  of  solar  declination 
and  Sun-Earth  distance  have  been  taken  for  the  middle  day  of 
each  of  the  months  January,  April,  July  and  October.  In 
all  cases  July  components  are  less  than  January  components  by 
about  7  per  cent  due  to  the  increase  in  Sun-Earth  distance, 
but  April  values  are  notably  different  from  solstitial 
values.  In  Fig.  4,  the  minima  in  April  components  at  45  km 
would  appear  to  arise  from  the  shape  of  the  latitudinal  profile 
of  the  Hough  1'unction  in  relation  to  that  of  the  heating  which 
decreases  towards  low  latitudes  as  discussed  above  (  §  6.1). 

The  differences  between  solstitial  and  April  values  in  Fig.  5 
are  again  most  likely  to  relate  to  changes  in  the  ozone 
distribution:  in  the  case  of  the  J^  _2  component  the  April 
values  are  actually  higher,  probably  as  a  consequence  cf  the 
form  of  which  unlike  that  for  other  modes  considered  has 

the  same  sign  at  all  latitudes.  The  results  plotted  in 
Figs.  4-  and  5  are  given  in  Table  5  together  with  the  October 


Height  profiles  of  Hough  components  of  heating 


Key: 


J1>n  for  n  =  1,3,5. 

- Jl,3; - (-)Ji,5’ 


(-)J 


1,1’ 


(-)J 


1,1 


from  Lindzen  (196?) 


values  which  hardly  differ  from  the  corresponding  April  values. 
Also  shown  on  Figs.  4  and  5  respectively  are  j}  ,  and  5 
taken  from  the  earlier  evaluation  of  Lindzen  (1967).  The 
new  values  are  considerably  lower  than  the  earlier  values  at 
heights  below  and  above  the  region  of  maximum  heating. 

The  largest  asymmetric  Hough  component  of  heating  is 
J^  ^  which  is  plotted  in  Fig.  6  for  January  and  October 
together  with  the  earlier  evaluation  of  Lindzen  (1967). 

Negative  values  arise  in  the  mesosphere  in  January  when 
N  hemisphere  ozone  densities  are  considerably  higher  than 
S  hemisphere  values  (Table  1)  leading  to  greater  heating 
in  spite  of  the  lower  wintertime  flux  of  incident  solar 
radiation.  In  the  stratosphere  in  January  the  asymmetry 
in  ozone  is  much  reduced  (Table  1)  and  the  heating  asymmetry 
conforms  with  that  of  the  solar  flux.  Likewise  the 
asymmetry  of  the  October  heating  accords  with  that  of  the 
radiational  flux,  the  solar  declination  being  -8°:  the 
October  asymmetry  in  the  stratosphere  increases  with  the 
seasonal  advance  until  the  January  profile  is  reached  while 
in  the  mesosphere  it  reverses.  July  and  April  profiles 
are  similar  to  those  shown  for  January  and  October  with  a 
reversal  of  sign.  The  higher  order  J?"  z,  J^  ^  components 
attain  maximum  values  of  17  and  9  mW/kg  respectively  and  are 
given  in  Table  3  whereas  the  asymmetric  components  J^  2 
J^  ^  are  less  than  2  and  1  mW/kg  respectively  and  are  not 
included. 
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mW/kg 


Fig.  6 


Height  profiles  of  the  component  of 

heating.  Key:  - present  calculations, 


N  hemisphere  winter  solstice  from 


Lindzen  (196?). 
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H 


(M.5) 


km 


mW/kq 


mW/kn 


Jan 

Apr 

Jul 

Oct 

Jan 

Apr 

3ul 

Oct 

n 

-0.0 

-0.0 

-0.0 

-0.0 

0.0 

-0.0 

-o.n 

0.0 

2 

-0.0 

-0.0 

-0.0 

-0.0 

0.0 

-0.0 

-0.0 

0.0 

4 

-0.0 

0.0 

-0.0 

-n.o 

0.0 

-o.o 

-0.0 

0.0 

6 

-0.0 

0.0 

-0.0 

0.0 

0.0 

-0.1 

-0.0 

0.1 

8 

0.0 

0.0 

0.0 

0.0 

0.1 

-0.1 

-0.1 

0.1 

10 

0.0 

0.0 

0.0 

0.0 

0.1 

-0.2 

-0.1 

0.2 

12 

0.0 

0.0 

0.0 

0.0 

0.1 

-0.3 

-0.1 

0.3 

14 

0.0 

0.0 

0.0 

0.0 

0.1 

-0.5 

-0.1 

0.4 

16 

0.0 

0.1 

0.0 

0.1 

0.3 

-0.7 

-0.3 

0.7 

18 

0.1 

0.1 

0.0 

0.1 

0.7 

-1.1 

-0.6 

1.0 

20 

-0.0 

-0.0 

-0.0 

-0.0 

1.5 

-1.5 

-1.4 

1.4 

22 

-0.1 

-0.2 

-0.1 

-0.2 

2.5 

-2.0 

-2.3 

2.0 

24 

-0.4 

-0.4 

-0.4 

-0.5 

3.8 

-2.6 

-3.4 

2.5 

26 

-0.9 

-0.8 

-0.8 

-0.8 

5.3 

-3.1 

-4.7 

3.1 

28 

-1.4 

-1.2 

-1.3 

-1.2 

7.1 

-3.8 

-6.5 

3.8 

30 

-2.0 

-1.5 

-1.9 

-1.6 

9.3 

-4.6 

-8.5 

4.6 

32 

-2.5 

-1.9 

-2.3 

-2.0 

12.5 

-6.4 

-11.5 

6.2 

34 

-3.1 

-2.4 

-2.9 

-2.5 

16.7 

-8.7 

-15.5 

a. 4 

36 

-3.7 

-3.0 

-3.4 

-3.0 

21.7 

-11.6 

-20.3 

11.0 

38 

-4.3 

-3.5 

-4.0 

-3.6 

27.0 

-14.6 

-25.4 

13.6 

40 

-5.1 

-4.2 

-4.7 

-4.3 

35.2 

-18.0 

-33.2 

16.4 

42 

-6.2 

-5.0 

-5.8 

-5.1 

44.4 

-21.9 

-42.1 

19.5 

44 

-6.8 

-4.7 

-6.4 

-4.7 

43.9 

-23.5 

-41.7 

20.4 

46 

-6.6 

-3.9 

-6.2 

-4.0 

33.7 

-19.4 

-32.2 

16.] 

48 

-5.8 

-3.6 

-5.5 

-3.6 

20.1 

-12.8 

-19.4 

9.P 

30 

-*.0 

-3.6 

-4.7 

-3.6 

6.5 

-7.6 

-6.8 

S.l 

52 

-4.3 

-3.7 

-4.0 

-3.7 

-5.1 

-5.6 

4.1 

3 .  R 

54 

-3.8 

-3.8 

-3.6 

-3.9 

-14.2 

-6.2 

12.7 

5.1 

56 

-3.4 

-4.0 

-3.2 

-4.0 

-20.8 

-8.7 

19.0 

8.1 

58 

-3.1 

-4.0 

-2.9 

-4.1 

-25.2 

-1°.Q 

23.1 

10.7 

60 

-2.7 

-3.9 

-2.6 

-3.9 

-27.0 

-12.1 

24.9 

12.0 

62 

-2.2 

-3.5 

-2.0 

-3.6 

-26.4 

-11.7 

24.4 

11.7 

64 

-1.6 

-2.9 

-1.5 

-3.0 

-23.5 

-10.3 

21.7 

10.3 

66 

-1.2 

-2.2 

-1.1 

-2.3 

-19.1 

-8.1 

17.7 

8.1 

68 

-0.8 

-1.5 

-0.8 

-1.6 

-14.3 

-5.8 

13.3 

5.7 

70 

-0.6 

-1.0 

-0.5 

-1.0 

-10.1 

-3.8 

9.4 

3.8 

72 

-0.4 

-0.6 

-0.4 

-0.6 

-6.6 

-2.2 

6.1 

2.1 

74 

-0.3 

-0.3 

-0.3 

-0.3 

-4.2 

-1.3 

3.8 

1.2 

76 

-0.3 

-0.2 

-0.2 

-0.2 

-2.5 

-0.8 

2.3 

0.7 

78 

-0.3 

-0.2 

-0.2 

-0.2 

-1.3 

-0.7 

1.1 

0.6 

80 

-0.3 

-0.4 

-0.3 

-0.4 

0.0 

-1.1 

-0.0 

1.1 

Table  3.  (continued) 


i 
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(1,1, -2) 

mW/kg 


(1,1, -3) 
mW/kq 


Jan 

Apr 

Jul 

Oct 

Jan 

Apr 

Jul 

Pet 

0 

-0.1 

-0.1 

-0.1 

-0.1 

-0.0 

n.n 

n.n 

-n.n 

2 

-0.1 

-0.1 

-0.1 

-P.l 

-0.0 

n.n 

n.n 

-n.n 

4 

-0.1 

-0.1 

-0.1 

-0.1 

-0.0 

n.n 

n.o 

-n.n 

6 

-0.2 

-0.2 

-0.2 

-0.2 

-0.0 

n.n 

n.n 

-n.n 

8 

-0.3 

-0.4 

-0.3 

-0.4 

-0.0 

n.i 

0.0 

-n.i 

10 

-0.5 

-0.6 

-0.5 

-0.6 

-0.1 

n.2 

n.i 

-n.i 

12 

-0.7 

-0.9 

-0.7 

-0.9 

-0.1 

0.2 

0.2 

-n.i 

14 

-1.0 

-1.3 

-1.0 

-1.3 

-0.3 

0.4 

0.3 

-n.2 

16 

-1.6 

-2.1 

-1.5 

-2.1 

-0.4 

0.5 

0.5 

-n.3 

18 

-2.9 

-3.7 

-2.7 

-3.7 

-0.6 

n.7 

0.7 

-0.3 

20 

-5.6 

-6.8 

-5.2 

-6.9 

-0.7 

0.9 

n.8 

-n.5 

22 

-10.0 

-11.4 

-9.3 

-11.6 

-0.6 

l.i 

n.8 

-n.6 

24 

-16.6 

-18.1 

-15.5 

-18.5 

-0.6 

1.4 

0.8 

-n.7 

26 

-25.0 

-26.0 

-23.2 

-26.6 

-0.6 

1.6 

0.9 

-i.n 

28 

-33.8 

-34.2 

-31.4 

-35.1 

-0.7 

1.9 

1.0 

-1.3 

30 

-42.6 

-42.6 

-39.7 

-43.8 

-1.0 

2.3 

J  .3 

-1.7 

32 

-53.4 

-53.7 

-49.8 

-55.1 

-1.3 

2.7 

1.7 

-2.2 

34 

-66.4 

-67.4 

-61.9 

-69.0 

-1.9 

3.3 

2.2 

-2.7 

36 

-B1.7 

-83.7 

-76.2 

-85.5 

-2.9 

4.1 

3.1 

-3.4 

38 

-99.7 

-103.0 

-93.0 

-105.1 

-4.4 

5.1 

4.5 

-4.3 

40 

-123.0 

-127.0 

-114.8 

-129.5 

-5.7 

6.1 

5.6 

-5.1 

42 

-151.7 

-156.8 

-141.5 

-159.7 

-7.0 

7.0 

6.9 

-5.8 

44 

-172.6 

-185.9 

-161.0 

-189.4 

-11.2 

10.4 

in. 6 

-8.6 

46 

-181.0 

-204.9 

-168.7 

-209.0 

-15.3 

14.7 

14.4 

-12.1 

48 

-181.1 

-211.6 

-168.6 

-216.1 

-17.5 

15. 

16.4 

-12.6 

SO 

-173.4 

-204.3 

-161.4 

-208.7 

-16.8 

13.7 

15.  R 

-in. 5 

52 

-162.6 

-186.9 

-151.2 

-190.9 

-14.2 

10.1 

13.5 

-7.3 

54 

-148.7 

-165.1 

-138.3 

-168.5 

-10.7 

8.2 

in. 3 

-6.5 

56 

-134.5 

-143.9 

-125.1 

-146.6 

-7.5 

8.6 

7.3 

-7.P 

58 

-119.2 

-125.0 

-110.9 

-127.2 

-4.8 

8.8 

4.8 

-8.6 

60 

-102.6 

-107.9 

-95.6 

-109.7 

-2.8 

7.9 

2.9 

-8.0 

62 

-83.1 

-91.7 

-77.4 

-93.1 

-1.4 

6.3 

1.6 

-6.3 

64 

-64.8 

-75.4 

-60.4 

-76.5 

-0.5 

4 .  R 

n.7 

-4.8 

66 

-48.2 

-58.8 

-44.9 

-59.7 

0.0 

3.7 

0.1 

-3.8 

68 

-34.8 

-43.5 

-32.4 

-44.2 

0.2 

3.1 

-0.1 

-3.1 

70 

-24.7 

-31.4 

-23.0 

-31.9 

0.3 

2.5 

-0.2 

-2.5 

72 

-17.3 

-21.4 

-16.1 

-21.8 

0.1 

1.8 

-0.0 

-1.8 

74 

-12.8 

-15.9 

-12.0 

-16.2 

-0.0 

1.3 

0.1 

-1.2 

76 

-10.3 

-13.1 

-9.6 

-13.4 

-0.1 

n.9 

n.i 

-0.8 

78 

-9.6 

-13.5 

-9.0 

-13.7 

-0.2 

0.8 

0.2 

-0.7 

80 

-11.1 

-17.9 

-10.3 

-18.3 

-0.3 

1.0 

0.3 

-0.9 

Table  3*  (continued) 


km 


mW/kg 

Jan  Apr  Jul  Oct 


mW/kg 

Jan  Apr  Jul  Oct 


n 

-0.0 

-0.0 

-0.0 

-0.0 

-0.0 

0.0 

0.0 

-n.o 

2 

-0.0 

-0.0 

-0.0 

-0.0 

-0.0 

0.0 

0.0 

-0.0 

4 

-0.0 

-0.0 

-0.0 

-0.0 

-0.0 

0.0 

0.0 

-0.0 

6 

-0.0 

0.0 

-0.0 

-0.0 

-0.0 

0.0 

0.0 

-o.n 

8 

-0.0 

0.0 

-0.0 

-o.n 

-0.0 

0.0 

0.0 

-n.o 

10 

-0.0 

n.i 

-0.0 

0.0 

-0.0 

n.o 

0.0 

-0.1 

12 

0.0 

n.2 

0.0 

0.1 

-0.1 

0.0 

0.0 

-0.1 

14 

0.2 

0.4 

0.2 

0.3 

-0.1 

-0.0 

0.1 

-0.1 

16 

0.3 

0.7 

0.3 

0.5 

-0.1 

-0.1 

0.1 

-n.o 

18 

0.2 

0.8 

0.3 

0.6 

-0.1 

-0.2 

0.1 

n.o 

20 

-0.3 

0.4 

-0.3 

0.1 

-0.1 

-0.2 

0.1 

0.1 

22 

-2.0 

-0.9 

-1.7 

-1.3 

-0.1 

-0.2 

0.1 

n.i 

24 

-4.9 

-3.3 

-4.3 

-3.8 

-n.i 

-0.1 

0.1 

-o.n 

26 

-8.6 

-6.4 

-7.8 

-6.9 

-0.2 

-0.0 

0.3 

-0.1 

28 

-12.3 

-9.3 

-11.2 

-9.8 

-0.4 

0.1 

0.4 

-0.2 

30 

-15.7 

-12.1 

-14.4 

-12.6 

-0.7 

0.1 

0.7 

-0.1 

32 

-19.6 

-15.8 

-18.0 

-16.3 

-0.9 

0.1 

0.9 

-0.1 

34 

-23.8 

-20.1 

-22.0 

-20.6 

-1.2 

0.1 

1.2 

-0.2 

36 

-28.3 

-24.6 

-26.2 

-25.3 

-1.6 

-0.0 

1.5 

-o.n 

38 

-33.4 

-29.2 

-30.9 

-29.8 

-1.8 

-0.5 

1.8 

0.4 

40 

-39.8 

-34.7 

-37.0 

-35.3 

-1.7 

-1.5 

1.7 

1.5 

42 

-48.4 

-41.8 

-45.0 

-42.4 

-1.7 

-2.3 

1.8 

2.3 

44 

-54.2 

-44.2 

-50.5 

-44.7 

-2.5 

-2.1 

2.5 

2.2 

46 

-36.4 

-38.6 

-52.8 

-38.5 

-3.2 

-A  .2 

3.3 

4.4 

48 

-56.4 

-29.7 

-53.0 

-28.8 

-3.5 

-8.4 

3.7 

R  .4 

30 

-54.2 

-22.1 

-51.1 

-20.0 

-3.4 

-9.S 

3.8 

8 .  R 

32 

-50.8 

-15.8 

-48.0 

-13.3. 

-2.5 

-7.6 

3.0 

6.6 

34 

-46.5 

-14.3 

-44.0 

-12.3 

-0.9 

-S.2 

1.4 

4.6 

36 

-42.5 

-17.«» 

-40.1 

-16.6 

1.4 

-3.7 

-0.8 

3.5 

38 

-37.6 

-22.3 

-35.5 

-21.6 

4.1 

-2.6 

-3.4 

2.6 

60 

-32.0 

-24.8 

-30.3 

-24.5 

6.1 

-1.7 

-5.4 

1.8 

62 

-25.2 

-24.5 

-23.8 

-24.5 

6.9 

-0.8 

-6.2 

0.9 

64 

-19.0 

-21.8 

-17.9 

-21.8 

6.7 

-0.3 

-6.0 

0.4 

66 

-13.5 

-16.9 

-12.8 

-17.0 

5.0 

-0.3 

-5.2 

0.3 

68 

-9.4 

-11.5 

-8.9 

-11.5 

4.5 

-0.5 

-4.1 

0.5 

70 

-6.5 

-6.9 

-6.2 

-6.9 

3.3 

-0.7 

-3.0 

0.7 

72 

-4.6 

-3.1 

-4.3 

-3.0 

2.2 

-0.8 

-2.0 

0.7 

74 

-3.4 

-1.0 

-3.3 

-O.o 

1.5 

-O.ft 

-1.3 

0.8 

76 

-2.8 

0.2 

-2.6 

0.4 

1.0 

-0.8 

-0.9 

0.7 

78 

-2.6 

0.4 

-2.4 

0.6 

0.7 

-0.8 

-0.7 

0.7 

80 

-2.9 

-0.6 

-2.8 

-0.4 

0.6 

-0.8 

-0.5 

0.8 

Table  5*  (continued) 
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H 

<2 

,2,5) 

1 7 

2.6) 

:* 

kn 

mW/kn 

mW/kcj 

Jan 

Apr 

Jul 

Hot 

Jan 

Apr 

Jul 

net 

n 

-o.n 

0.0 

n.n 

0.0 

0.0 

0.0 

0.0 

o.n  i 

2 

-n.o 

n.n 

o.n 

n.o 

0.0 

0.0 

0.0 

n.n  .i 

4 

-n.n 

n.n 

0.0 

n.n 

0.0 

0.0 

n.n 

n.o  ) 

6 

-0.0 

0.0 

n.o 

n.o 

0.0 

0.0 

0.0 

n.n 

R 

-0.0 

n.n 

0.0 

n.o 

n.n 

n.o 

o.n 

n.n 

in 

-n.i 

0.0 

n.i 

n.o 

0.1 

0.1 

n.o 

n.i 

i? 

-0.1 

n.n 

n.i 

n.n 

n.i 

0.1 

n.i 

n.i 

14 

-0.2 

o.o 

0.2 

-0,0 

0.1 

0.1 

0.1 

n.2  i'i 

16 

-0.2 

0.0 

0.2 

-n.o 

0.2 

0.2 

0.2 

0.2  •, 

1ft 

-0.4 

0.1 

0.3 

-n.o 

0.3 

0.4 

0.2 

0.4  K 

20 

-0.5 

0.1 

0.5 

-n.i 

0.4 

0.5 

0.3 

0.6 

22 

-o.o 

0.2 

0.8 

-0,1 

0.5 

0.8 

0.4 

0.8  !1 

24 

-1.3 

0.4 

1.2 

-0.2 

0.7 

l.n 

0.6 

l.n  ?! 

26 

-1.6 

0.5 

1.6 

-0.3 

0.9 

1.4 

0.8 

1.4  i- 

28 

-1.8 

0.7 

1 .8 

-0.5 

1.1 

1.8 

0.9 

1.8 

30 

-1.8 

0.9 

1.8 

-0.6 

1.2 

2.3 

1.0 

2.3 

32 

-1.9 

1.0 

1.9 

-0.7 

1.5 

3.0 

1.2 

3.0  *  : 

34 

-2.0 

1.1 

2.0 

-0.7 

1.9 

3.9 

1.7 

36 

-2.2 

1.0 

2.1 

-0.7 

2.5 

5.1 

2.2 

5.2  ‘  1 

38 

-2.7 

1.2 

2.7 

-0.9 

3.3 

6.8 

2.9 

7.1 

40 

-3.2 

1.6 

3.1 

-1.3 

4.7 

8.7 

4.2 

9.0 

42 

-2.9 

1.2 

2.8 

-l.n 

6.3 

]  I  .2 

5.6 

11.5 

44 

-5.5 

2.4 

5.2 

-2.1 

6.3 

15.0 

5.6 

16.5 

46 

-11.3 

6.4 

in. 6 

-5.9 

6.6 

1°.R 

5.9 

20.7 

4R 

-18.0 

10.1 

16.0 

-9.2 

7.6 

10.4 

6.7 

20 . 3 

5n 

-21.o 

11.2 

20.6 

-m.n 

8.4 

15.7 

7.4 

16.2 

52 

-23.6 

in.  s 

22.  1 

-o.? 

8.5 

IK! 

7.6 

11.5 

54 

-23.3 

o.O 

21 

-7.o 

8.5 

7.5 

7.7 

7.7 

56 

-21.8 

7.4 

20.5 

-6.4 

R.o 

5.0 

8.1 

5.2 

58 

-10.6 

5.8 

18.4 

-4.9 

0.6 

3.6 

9.8 

3.8  I 

60 

-17.0 

4.3 

16.0 

-3.7 

9.0 

2.o 

9.2 

3.1 

6. 

-14.2 

3.2 

13.3 

-2.7 

9.6 

2.6 

8.9 

2.8  v 

64 

-11.4 

2.4 

10.8 

-1.9 

8.5 

2.5 

7.9 

2.6 

66 

-B.B 

1.8 

8.3 

-1.4 

7.0 

2.3 

6.5 

2.4 

68 

-6.5 

1.3 

6.1 

-1.1 

•  5.3 

2.0 

4.9 

2.1 

70 

-4.6 

1.0 

4.4 

-0.8 

3.9 

1.8 

3.6 

1.9 

72 

-3.2 

0.7 

3.0 

-0.6 

2.7 

1.6 

2.5 

1.6  |! 

74 

-2.2 

0.6 

2.1 

-0.5 

1.9 

1.5 

1.8 

K5 

76 

-1.6 

0.5 

1.5 

-0.5 

1.4 

1.4 

1.3 

1.4 

78 

-1.2 

0.6 

1.2 

-0.5 

1.1 

1.4 

1.0 

1.4 

80 

-1.3 

0.8 

1.2 

-0.8 

0.9 

1.5 

0.8 

1.5 

Table  3.  (continued) 


The  results  obtained  for  symmetric  Hough  components  of 
2  2  2 

heating;  J2,2’  J2  4’  J2  6  are  shown  in  FiS*  7.  Semi-diurnal 
Hough  functions  for  s  =  2  have  previously  been  presented 
graphically  (Groves,  1975)  but  may  differ  in  sign  from  ti.ose 
used  here  which  are  positive  at  the  equator  if  equatorially 
symmetric  or  increase  at  the  equator  with  latitude  if 
equatorially  asymmetric.  As  with  diurnal  components  the 

imaginary  oarts  of  J~  are  relatively  small  (of  order  10  7) 

»n 

and  will  therefore  be  neglected.  July  values  are  slightly 

less  than  January  values  due  to  the  increase  of  the  Sun-3arth 

distance:  October  values  are  not  shown  being  nearly  the 

p 

same  as  April  values.  J0  ■>  shows  very  little  annual 
variation  due  possibly  to  the  form  of  0,  „  which  unlike 
that  for  other  modes  considered  has  the  same  sign  at  ail 
latitudes.  Fig.  7  shows  that  at  many  heights  the  present 
values  of  J ^  ^  are  considerably  lower  tr.ar.  ti.ose  of  Chapman 
and  Lindzen  (1970)  and  that  they  compare  puite  well  with  those 

P  " 

of  Butler  and  Small  (1965)*  J-,  have  April  and 

October  values  greatly  exceeding  January  and  July  values  in 
the  region  of  45  km  whereas  in  the  mesosphere  January 
and  July  values  greatly  exceed  April  and  October  values: 
such  differences  are  intricately  related  to  the  form  of  the 
ozone  distribution  and  that  of  the  Hough  functions. 


f 
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Fig.  8  shows  the  asymmetric  componer.ts  J 


2 

2,3’ 


J 


2 

2,5 


calculated  for  January  and  October;  the  July  and  April 
profiles  differing  only  slightly  from  these  with  signs 
reversed.  At  the  solstices,  semi-diurnal  heating  is 
greater  in  the  winter  than  the  summer  hemisphere  as  the 
shorter  duration  of  daylight  generates  a  larger  Fourier 
component.  This  difference  is  further  enhanced  above  45  km 
by  a  large  seasonal  asymmetry  in  ozone  densities  (Table  1); 
and  the  rapid  increase  of  beating  rates  with  height  apparent 
in  Fig.  8  results.  Semi-diurnal  components  are  tabulated 
in  Table  3* 


7.  Discussion 

Formulae  have  been  derived  in  %3  for  the  Hough  components 
of  ozone  heating  in  the  presence  of  a  lower  reflecting  layer, 
account  being  taken  in  the  theory  of  longitudinal  dependences 
in  the  ozone  distribution  and  the  reflective  properties  of 
the  underlying  layer.  A  simple  representation  of  surface 
and  cloud  reflectivities  has  been  introduced  ( ^  5)  which 
depends  on  longitude,  but  insufficient  data  are  available 
for  incorporating  a  longitudinal  dependence  into  the  ozone 
model  introduced  in  §4.  Accordingly  the  Hough  components 
presented  here  are  those  of  migrating  modes,  i.e.  those  for 
which  s  *  m,  as  these  depend  only  on  mean  longitudinal 
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properties . 

Contributions  from  longitudinal  variations  in  lower 
layer  reflectivity,  mainly  in  the  cloud  distribution,  have 
been  investigated  as  total  cloudiness  increases  stratospheric 
heating  by  as  much  as  20  per  cent  at  30  km  (§  6.1);  but  the 
Hough  components  obtained  for  s  ^  m  were  mostly  much  less 
than  3  per  cent  of  the  corresponding  s  *  m  components  at  the 
same  height  and  have  therefore,  not  been  presented. 

Furthermore  such  components  maximize  at  about  32  km  and 
would  be  significantly  affected  by  any  longitudinal  ozone 
variations  of  the  magnitude  (12  per  cent)  previously 
indicated  (5  4) . 

A  latitudinal  and  seasonal  model  of  ozone  densities 
(Table  1)  has  been  adopted  which  is  based  on  observational 
data  in  contrast  to  the  1964  photochemical  model  of  Leovy 
which  was  used  in  earlier  Hough  component  calculations. 
Although  the  available  ozone  data  have  obvious  limitations, 
particularly  at  high  latitudes  and  above  50  km,  the  model 
is  expected  to  be  a  useful  improvement  on  the  1964  theoretical 
model  which  has  its  own  limitations,  such  as  that  of  being 
based  on  pure  oxygen  reactions.  For  a  hydrogen-oxygen- 
nitrogen  system  lower  ozone  concentrations  are  computed  in 
the  mesosphere  due  to  hydrogen-oxygen  reactions  and  also 
in  the  stratosphere  due  to  nitrogen-oxygen  reactions 
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(Shimazaki  and  Wuebbles,  1973;  Park  and  London,  1974).  In 
the  lower  stratosphere  departures  from  radiational- 
photochemical  equilibrium  occur  as  ozone  lifetimes  increase 
sufficiently  for  vertical  mixing  to  become  an  effective 
removal  process,  reducing  concentrations  still  further. 

Comparisons  have  been  made  in  Figs.  4,  3  and  7  respectively 
between  the  newly  calculated  Hough  components  _p 

p 

and  J2  2  and  tllose  of  earlier  calculations  (Lindzen,  1967; 
Chapman  and  Lindzen,  1970):  the  new  values  are  considerably 
letter  than  the  earlier  values  in  height  regions  above  and 
below  the  stratopause  region  of  maximum  heating  in  reasonable 
accord  with  the  introduction  of  lower  o zone  densities.*  The 

p 

newly  calculated  J2  p  have  been  compared  in  Fig.  7  with  the 
evaluation  of  But] er  and  Small  (1963)  which  was  based  on  a 
simple  ozone  model  devised  from  observational  data:  in  view 
of  the  limited  ozone  and  solar  spectral  UV  data  available  at 
that  time,  the  agreement  between  the  two  sets  of  results  is 
better  than  expected. 

In  general  the  value  of  any  Hough  component  depends 
intricately  on  the  relation  between  the  latitudinal  profile 
of  ozone  and  that  of  the  particular  Hough  function.  The 
decrease  in  April  ozone  densities  towards  low  latitudes 
between  40  and  50  km  (Table  1)  would  appear  to  be  responsible 
for  the  minimum  at  48  km  in  the  components  in  Fig.  4;  and 

*  Since  this  paper  was  finalized,  reference  has  been  obt  ned  to  Forbes, 
J.M.  and  Garrett,  H.B.  (1978)  Geophys.  Res.  Lett.  £,  1013-  ul6  reporting 
similar  significant  differences  from  the  heating  profiles  given  by  Chapman 
and  Lindzen  (1970). 
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2  2 

also  for  the  higher  April  values  of  ^2  6  ^S*  5» 

the  effect  being  in  the  opposite  sense  due  to  the  different 
form  of  the  Hough  functions  in  the  two  cases. 

With  regard  to  asymmetric  Hough  components  a  large  part 
of  the  solstitial  asymmetry  in  the  mesosphere  (Figs.  6  and  8) 
is  generated  by  the  asymmetry  in  mesospheric  ozone  densities 
shown  in  Table  1  which  has  winter  values  greatly  in  excess 
of  summer  values,  being  in  the  opposite  sense  to  that  predicted 
by  photochemical  theory  (Park  and  London,  1974).  This 
disagreement  between  observation  and  theory  has  previously 
been  reported  by  Watanabe  and  Tohmatsu  (1976)  on  the  evidence 
of  a  series  of  rocket-borne  solar  UV  radiometer  measurements 
conducted  since  1965  at  31°  15 'N,  150°  05' E;  a  maximum 
winter-summer  ratio  of  2.8  being  obtained  at  65  km  which 
is  comparable  with  that  shown  by  Table  1.  Observational 
data  above  50  km  are  sparse  with  values  at  the  same  latitude 
sometimes  differing  by  a  factor  of  2  or  more  and  hence  the 
results  in  this  paper  above  50  km  are  particularly  open  to 
review  in  the  light  of  further  data  on  ozone  densities. 
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